LINER/H ii "Transition" Nuclei and the Nature of NGC 4569 
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ABSTRACT 

Motivated by the discovery of young, massive stars in the nuclei of some LINER/H II "transi- 
tion" nuclei such as NGC 4569, we have computed photoionization models to determine whether 
some of these objects may be powered solely by young star clusters rather than by accretion- 
powered active nuclei. The models were calculated with the photoionization code CLOUDY, 
using evolving starburst continua generated by the the STARBURST99 code of Leitherer et 
al. (1999). We find that the models are able to reproduce the emission-line spectra of transi- 
tion nuclei, but only for instantaneous bursts of solar or higher metallicity, and only for ages of 
^3 — 5 Myr, the period when the extreme-ultraviolet continuum is dominated by emission from 
Wolf-Rayet stars. For clusters younger than 3 Myr or older than 6 Myr, and for models with a 
constant star-formation rate, the softer ionizing continuum results in an emission spectrum more 
typical of H II regions. This model predicts that Wolf-Rayet emission features should appear in 
the spectra of transition nuclei. While such features have not generally been detected to date, 
they could be revealed in observations having higher spatial resolution. Demographic arguments 
suggest that this starburst model may not apply to the majority of transition nuclei, particularly 
those in early- type host galaxies, but it could account for some members of the transition class 
in hosts of type Sa and later. The starburst models during the Wolf-Rayet-dominated phase can 
also reproduce the narrow-line spectra of some LINERs, but only under conditions of above-solar 
metallicity and only if high-density gas is present (n e > 10 5 cm~ 3 ). This scenario could be 
applicable to some "Type 2" LINERs which do not show any clear signs of nonstellar activity. 
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Introduction 



Emission-line nebulae in galactic nuclei are gen- 
erally considered to fall into three major cate- 
gories: star-forming or H II nuclei, Seyfert nu- 
clei, and low-ionization nuclear emission-line re- 
gions, or LINERs. The formal divisions between 
these classes are somewhat arbitrary, as the ob- 
served emission-line ratios of nearby galactic nu- 
clei fall in a continuous distribution between LIN- 
ERs and Seyfert nuclei and between LINERs and 
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H II nuclei (e.g., Ho et al. 1993). Traditionally, 
LINERs have been defined as those nuclei hav- 
ing emission-line flux ratios which satisfy the re- 
lations [O II] A3727/[0 III] A5007 > 1 and [O I] 
A6300/[O III] A5007 > 1/3 (Heckman 1980). It 
is possible to construct alternative but practically 
equivalent definitions, based on other line ratios, 
which can be applied to datasets that do not in- 
clude the wavelengths of the [O I], [O II], or [O III] 
lines (e.g., Ho et al. 1997a). 
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A sizeable minority of galactic nuclei has 
emission-line ratios which are intermediate be- 
tween those of "pure" LINERs and those of typical 
H II regions powered by hot stars; these galaxies 
would be classified as LINERs except that their 
[O I] A6300 line strengths are too small in compar- 
ison with other lines to meet the formal LINER 
criteria. Objects falling into this category have 
been dubbed "transition" galaxies by Ho et al. 
(1993), and although this nomenclature is some- 
what ambiguous we adopt it here for consistency 
with the spectroscopic survey of Ho et al. (1997a). 
That survey defined the transition class in terms 
of the following flux ratios: 

[O III] A5007/H/3 < 3, 
0.08 < [O I] A6300/Ha < 0.17, 

[N II] A6583/Ha > 0.6, 
[S II] AA6716, 6731 /Ha > 0.4. 

Filippcnko & Tcrlevich (1992) have used the term 
"weak-[0 I] LINERs" to refer to galaxies having 
[N II] A6583/Ha > 0.6 (typical of LINERs) but 
which have [O I] A6300/Ha < 1/6. This cate- 
gory is essentially identical to the transition class 
of Ho et al. (1993, 1997a), and we will refer to 
these galaxies as transition objects in this paper. 

According to the survey results of Ho et al. 
(1997b), this transition class accounts for 13% of 
all nearby galaxies, making them about as numer- 
ous as Seyfert nuclei. The Hubble type distri- 
bution of transition galaxies is intermediate be- 
tween that of LINERs, which are most common in 
E/SO/Sa galaxies, and that of H II nuclei, which 
occur most often in Hubble types later than Sb 
(Ho et al. 1997b). Roughly 20% of galaxies with 
Hubble types ranging from SO to Sbc belong to 
the transition class. There is not a consensus, how- 
ever, as to whether these transition objects should 
be regarded as star-forming nuclei, as accretion- 
powered active nuclei, or as composite objects 
powered by an AGN and by hot stars in roughly 
equal proportion. 

There is a large body of literature on the sub- 
ject of the excitation mechanism of LINERs which 
is relevant to the similar transition class. A vari- 
ety of physical mechanisms has been proposed to 
explain the emission spectra of LINERs, including 
shocks, photoionization by a nonstellar ultraviolet 
(UV) and X-ray continuum, and photoionization 
by hot stars. (See Filippenko 1996 for a review.) 



The possibility that LINERs (and Seyfert nuclei as 
well) might be photoionized by starlight was raised 
by Terlevich & Melnick (1985), who suggested that 
very hot (T off - 10 5 K) Wolf-Rayet (W-R) stars 
in a metal-rich starburst could give rise to an ion- 
izing continuum with a nearly power-law shape in 
the extreme-UV. More recent atmosphere mod- 
els have indicated substantially lower tempera- 
tures for W-R stars, however, casting doubt on the 
Warmer hypothesis (Leitherer et al. 1992). Subse- 
quent photoionization models have attempted to 
explain LINER and transition-type spectra as re- 
sulting from massive main-sequence stars. Filip- 
penko & Terlevich (1992) found that the spectra 
of weak-[0 I] LINERs could be explained in terms 
of photoionization by 03-04 stars having effec- 
tive temperatures of > 45, 000 K, at ionization pa- 
rameters of U « 10~ 3 - 7 to 10- 3 - 3 . Shields (1992) 
carried this line of argument farther, proposing 
that genuine LINER spectra could be generated by 
early O stars with T e g ~ 50, 000 K, provided that 
a high-density component (n e « 10 5 5 cm~ 3 ) is 
present in the NLR; the high densities are needed 
to boost the strengths of high critical-density emis- 
sion lines, most notably [O I] A6300. Similar con- 
clusions were reached by Schulz & Fritsch (1994), 
who explored the effects of absorption by ionized 
gas as a means to harden the effective ionizing 
spectrum. Recent observations, particularly in the 
UV and X-ray bands, have provided convincing 
evidence that many LINERs are in fact AGNs, 
particularly the "Type 1" LINERs which have a 
broad component to the Ha emission line (for a 
recent review see Ho 1999). The possibility has 
remained, however, that some LINERs and tran- 
sition nuclei are powered entirely by bursts of star 
formation. 

An important shortcoming of the model cal- 
culations performed by Filippenko & Terlevich 
(1992) and Shields (1992) is that the ionizing con- 
tinua used as input were those of single O-type 
stars; these studies did not address the question 
of whether a LINER or transition-type spectrum 
could result from the the integrated ionizing con- 
tinuum of a young stellar cluster. Compared with 
these single-star models, the contribution of late- 
O and B stars will soften the ionizing spectrum, 
making the emission-line ratios tend toward those 
of normal H II regions. W-R stars, on the other 
hand, will harden the ionizing spectrum during 



2 



the period when these stars are present, roughly 
3 — 6 Myr after the burst. Another drawback of 
the O-star models is that they require the pres- 
ence of stars with effective temperatures higher 
than are thought to occur in H II regions of so- 
lar or above-solar metallicity, in order to produce 
a LINER or transition-type spectrum rather than 
an H II region spectrum. Their applicability to 
galactic nuclei is therefore somewhat unclear. 

Other mechanisms have been proposed for gen- 
erating LINER or transition-type spectra. Shock 
excitation by supernova remnants in an aging star- 
burst may give rise to some transition objects; the 
nucleus of NGC 253 is a likely candidate for such 
an object (Engelbracht et al. 1998). Also, post- 
AGB stars and planetary nebula nuclei will pro- 
duce a diffuse ionizing radiation field which could 
be responsible for the very faint LINER emission 
(with Ha equivalent widths of ~ 1 A) observed in 
some ellipticals and spiral bulges (Binette et al. 
1994). 

An alternate possibility is that the transition 
galaxies may simply be composite systems con- 
sisting of an active nucleus surrounded by star- 
forming regions. For a galaxy at a distance of 10 
Mpc, for example, a 2"-wide spectroscopic aper- 
ture will include H II regions within 50 pc of 
the nucleus. Galaxies having emission lines both 
from a LINER nucleus and from surrounding star- 
forming regions, in roughly equal proportions, will 
appear to have a transition-type spectrum. This 
interpretation was advocated by Ho et al. (1993) 
as the most likely explanation for the majority of 
transition galaxies, and is consistent with the ob- 
served Hubble type distribution for the transition 
class. Other authors have similarly contended that 
transition galaxies are AGN/H II region compos- 
ites, based on optical line-profile decompositions 
(Veron et al. 1997; Goncalves et al. 1999) and 
near-infrared spectra (Hill et al. 1999). Two of 
the 65 transition nuclei observed in the Ho et al. 
(1997b) survey have a broad component to the Ha 
emission line, indicating the likely presence of an 
AGN, and it is probable that many more tran- 
sition nuclei contain obscured AGNs which were 
not detected in the optical spectra. On the other 
hand, radio observations do not appear to support 
the composite AGN/starburst interpretation. In a 
VLA survey of nearby galactic nuclei, Nagar et al. 
(1999) find compact, flat-spectrum radio cores in 



more than 50% of LINER nuclei, but in only 6% 
(1 of 18) of transition objects. This discrepancy 
suggests that the simple picture of an ordinary 
LINER surrounded by star-forming regions may 
not apply to the majority of transition objects. 

Recent results from the Hubble Space Telescope 
(HST) have shed new light on the question of 
the excitation mechanism of transition nuclei. As 
shown by Maoz et al. (1998), the UV spectrum 
of the well-known transition nucleus in NGC 4569 
over 1200-1600 A is virtually identical to that of 
a W-R knot in the starburst galaxy NGC 1741, 
indicating that O stars with ages of a few Myr 
dominate the UV continuum. Maoz et al. (1998) 
find that the nuclear star cluster in NGC 4569 
is producing sufficient UV photons to ionize the 
surrounding narrow-line region, a key conclusion 
which provides fresh motivation to study stellar 
photoionization models. The brightness of the 
NGC 4569 nucleus, and the consequently high S/N 
observations that have been obtained, make it one 
of the best objects with which to study the tran- 
sition phenomenon. 

The recent availability of the STARBURST99 
model set (Leitherer et al. 1999) has prompted us 
to reexamine the issue of ionization by hot stars in 
LINERs and transition nuclei. These models give 
predictions for the spectrum and luminosity of a 
young star cluster, for a range of values of clus- 
ter age, metal abundance, and stellar initial mass 
function (IMF) properties. Using the photoioniza- 
tion code CLOUDY (Ferland et al. 1998) in com- 
bination with the STARBURST99 model con- 
tinua, we have calculated the expected emission- 
line spectrum of an H II region illuminated by 
a young star cluster, to test the hypothesis that 
some LINERs and transition nuclei may be pow- 
ered by starlight. Similar calculations have been 
performed by Stasinska & Leitherer (1996), but 
for the physically distinct case of metal-poor ob- 
jects representing H II galaxies. Other examples 
of photoionization calculations for H II regions us- 
ing evolving starburst continua are presented by 
Garcia- Vargas & Diaz (1994), Garcia- Vargas et al. 
(1995), and Bresolin et al. (1999). 

2. The Nucleus of NGC 4569 

Before describing the photoionization model- 
ing, we review the properties of NGC 4569, as it is 
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among the best-known examples of the transition 
class. NGC 4569 is a Virgo cluster spiral of type 
Sab, with a heliocentric velocity of —235 km s _1 , 
and we assume a distance of 16.8 Mpc for consis- 
tency with the catalog of Ho ct al. (1997a). Its nu- 
cleus is remarkably bright for a non-Seyfert, and 
so compact in the optical that Humason (1936) 
suspected it to be a foreground Galactic star. It is 
also an unusually bright UV source, with the high- 
est 2200 A luminosity of the LINERs and transi- 
tion objects observed by Maoz et al. (1995) and 
Barth et al. (1998). HST Faint Object Spectro- 
graph (FOS) spectra show that the UV continuum 
is dominated by massive stars, with prominent P 
Cygni profiles of C IV A1549, Si IV A1400, and 
N V A1240 (Maoz et al. 1998). The UV spectrum 
is nearly an exact match to the spectrum of one of 
the starburst knots in the W-R galaxy NGC 1741, 
an object with a likely age in the range 3-6 Myr 
(Conti et al. 1996). The optical spectrum of the 
NGC 4569 nucleus is dominated by the light of 
A-type supergiants, providing additional evidence 
for recent star formation (Keel 1996). 

One key result of the Maoz et al. (1998) study 
was the conclusion that the nuclear starburst in 
NGC 4569 is producing sufficient numbers of ion- 
izing photons to power the narrow-line region, as- 
suming that the surrounding nebula is ionization- 
bounded, even without correcting for the effects 
of internal extinction on the UV continuum flux. 
In fact, there appears to be substantial extinction 
within NGC 4569, as demonstrated by the UV 
continuum slope as well as the presence of deep 
interstellar absorption features (Maoz et al. 1998). 
Ho et al. (1997a) derive an internal reddening of 
E(B-V) = 0.46 mag from the Ha/H/3 ratio, while 
Maoz ct al. (1998) estimate a UV extinction of 
A w 4.8 mag at 1300 A by comparison of the ob- 
served UV slope with the expected spectral shape 
of an unreddened starburst. 

Despite the fact that NGC 4569 is often referred 
to as a LINER, and in some cases presumed to 
contain an AGN on the basis of that classification, 
there is no single piece of evidence which conclu- 
sively demonstrates that an AGN is in fact present 
at all. The HST images and spectra are all con- 
sistent with the nucleus being a young, luminous, 
and compact starburst region. No broad-line com- 
ponent is detected on the Ha emission line (Ho 
et al. 1997a), and no narrow or broad emission 



lines are visible at all in the UV spectrum other 
than the P Cygni features that are generated in 
O-star winds (Maoz et al. 1998). Only the optical 
emission-line ratios point to a possible AGN clas- 
sification. In a thorough study of optical and IUE 
UV spectra, Keel (1996) concluded that there was 
at best weak evidence for the presence of an AGN 
in NGC 4569, and that any AGN continuum com- 
ponent, if present, must have an unusually steep 
spectrum. 

Furthermore, while the nucleus of NGC 4569 is 
certainly extremely compact, the UV and optical 
HST images show that the nucleus is not domi- 
nated by a central point source. At 2200 A, the 
nucleus appears extended in WFPC2 images with 
FWHM sizes of 13 and 9 pc along its major and 
minor axes (Barth et al. 1998). Optical WFPC2 
images have been discussed recently by Pogge et 
al. (1999), who state that the nucleus is unresolved 
by HST. We have obtained these same images from 
the HST archive. While the nucleus is certainly 
compact, we find that it is clearly extended even at 
the smallest radii. A 12-second, CR-SPLIT expo- 
sure in the F547M (U-band) filter is unsaturated 
and allows a radial profile measurement. We find 
a FWHM size of 14 pc by 8 pc along the major and 
minor axes, consistent with the size of the nuclear 
cluster measured at 2200 A. Barth et al. (1998) es- 
timated that at most 23% of the nuclear UV flux 
could come from a central point source. From the 
equivalent widths of stellar-wind features in the 
UV spectrum, Maoz et al. (1998) give a similar 
upper limit of ~ 20% to the possible contribution 
of a truly featureless continuum to the observed 
UV flux. 

X-ray observations of NGC 4569 with ROSAT 
have revealed a source coincident with the nucleus 
which is unresolved at the 2" resolution of the HRI 
camera (Colbert & Mushotzky 1999). This does 
not necessarily indicate that an AGN is present, 
however, as the optical/UV size of the starburst 
core is an order of magnitude smaller than the 
HRI resolution. ASCA observations show that the 
X-ray emission is extended over arcminute scales 
in both the hard (2-7 keV) and soft (0.5-2 keV) 
bands (Terashima et al. 1999). Interestingly, the 
compact source seen in the ROSAT image is de- 
tected only in the soft ASCA band, while there is 
no detectable contribution from a compact, hard 
X-ray source. The spectral shape of the compact 
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soft X-ray component is consistent with an origin 
either in an AGN or in X-ray binaries (Tschokc 
& Hensler 1999), but the lack of a compact hard 
X-ray source argues against the AGN interpreta- 
tion. If an AGN is present, it must be highly ob- 
scured even at hard X-ray energies, with an ob- 
scuring column of Nh > 10 23 cm~ 2 (Tcrashima 
et al. 1999). In radio emission, VLA observations 
show that the NGC 4569 nucleus is an extended 
source with a size of 4" and no apparent core (Neff 
& Hutchings 1992), in contrast with the compact, 
AGN-like cores found in some LINERs (Falcke et 
al. 1998). 

Shock excitation has often been considered as 
a mechanism to power the narrow emission lines 
in LINERs. However, the lack of narrow emis- 
sion features in the UV spectrum of NGC 4569 
argues against shock-heating models for this ob- 
ject, as existing shock models generally predict 
strong UV line emission (e.g., Dopita & Suther- 
land 1996). Shock-excited filaments in supernova 
remnants show strong emission in high-excitation 
UV lines such as C IV A1549 and He II A1640 
(e.g., Blair et al. 1991, 1995) which are altogether 
absent from the NGC 4569 spectrum. Similarly, 
the shock-excited nuclear disk of M87 (Dopita et 
al. 1997) has a high-excitation UV line spectrum 
which bears no resemblance to the NGC 4569 
spectrum. From an analysis of infrared spectra, 
Alonso-Herrero et al. (1999) proposed that the 
NGC 4569 nucleus is powered by an 8-11 Myr-old 
star burst, by a combination of stellar photoioniza- 
tion and shock heating from supernova remnants. 
While this hypothesis may be applicable to some 
LINERs and transition galaxies, the UV spectrum 
of NGC 4569 shown by Maoz et al. (1998) is incon- 
sistent with a burst of such advanced age, as the 
P Cygni features of C IV, Si IV, and N V would 
have disappeared from a single-burst population 
after about 6 Myr. 

The overall picture emerging from these obser- 
vations is that the NGC 4569 nucleus is a com- 
pact, luminous, and young starburst. The only 
reason to invoke the presence of an AGN at all 
would be to explain the higher strengths of the 
low-ionization forbidden lines in comparison with 
values observed in normal H II nuclei. If it were 
indeed possible for a young starburst to produce 
transition or LINER-type emission lines in the sur- 
rounding gas, then there would be no reason to 



consider AGN models for NGC 4569. 

3. Photoionization Calculations 
3.1. The Ionizing Continuum 

As discussed by Filippenko & Terlevich (1992) 
and Shields (1992), the key ingredient neces- 
sary for generating a LINER or transition-type 
emission-line spectrum is an ionizing continuum 
which is harder than that produced by typical 
clusters of OB stars. A harder continuum will 
produce a more extended partially-ionized zone in 
the surrounding H II region, boosting the strength 
of the low-ionization lines which are typical of 
LINER spectra: [O I] A6300, [O II] A3727, [N II] 
AA6548,6583, and [S II] AA6716,6731. 

To represent the ionizing continuum of a young 
starburst, we have chosen the STARBURST99 
model set; we refer the reader to Lcitherer et 
al. (1999) for the details of the methods used 
to construct these models. Briefly, the STAR- 
BURST99 code employs the Geneva stellar evolu- 
tion models of Meynet et al. (1994), with enhanced 
mass-loss rates, for high-mass stars. Atmospheres 
are represented by the models compiled by Lejeune 
et al. (1997) and Schmutz et al. (1992). Figures 1- 
12 of Leitherer et al. (1999) display the spectral en- 
ergy distributions of the STARBURST99 model 
clusters for a range of burst ages and for a variety 
of initial conditions. From the figures, some im- 
portant trends are readily apparent. During the 
first 2 Myr after an instantaneous burst, the con- 
tinuum is dominated by the hottest O stars, and 
there is essentially no emission below 228 A, cor- 
responding to the ionization energy of Hc + . The 
appearance of W-R stars during the period 3-5 
Myr after the burst results in a dramatic change 
in the UV continuum, as these stars emit strongly 
in the He ++ continuum below 228 A. From 6 Myr 
onwards, the W-R stars disappear and the UV 
continuum rapidly fades and softens as the burst 
ages. Only the models with an upper mass limit of 
-^up = 100M Q generate the hard, W-R-dominated 
UV continuum; the model sequences with M up = 
30M Q do not generate significant numbers of pho- 
tons below 228 A for any ages because the pro- 
genitors of W-R stars are not present in the initial 
burst. Constant star- formation rate models with 
Map = 100M Q form W-R stars continuously after 
3 Myr, but the overall shape of the UV continuum 
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is softer than in the instantaneous burst models, 
because of the continuous formation of luminous 
O stars. 

These results provide a useful starting point for 
the photoionization calculations. If it is possible 
for the H II region surrounding a young cluster 
to resemble a LINER or transition object, then 
this is most likely to occur when the ionizing con- 
tinuum is hardest, when W-R stars are present 
during t « 3 — 5 Myr after a burst. Very massive 
stars (in the range 30-100 M Q or greater) must 
be present in the burst or else the requisite W- 
R stars will not appear. The formation of W-R 
stars is enhanced at high metallicity, so the ability 
to generate a LINER or transition-type spectrum 
may be a strong function of metal abundance as 
well as age. 

3.2. Model Grid 

To create grids of photoionization models, we 
fed the UV continua generated by the STAR- 
BURST99 models into the photoionization code 
CLOUDY (version 90.04; Ferland et al. 1998). For 
each time step, a grid of models was calculated by 
varying the nebular density and the ionization pa- 
rameter, which is defined as the ratio of ionizing 
photon density to the gas density at the ionized 
face of a cloud. Real LINERs and transition nu- 
clei are likely to contain clouds with a range of 
values of density and ionization parameter, and 
more general models incorporating density and 
ionization stratification can be constructed as lin- 
ear combinations of these simple single-zone mod- 
els. 

All models were run with the following range 
of parameters: burst age from 1 to 10 Myr at in- 
crements of 1 Myr, with log U ranging from —2 
to —4 at increments of 0.5, and a constant density 
ranging from log (nn/cm" 3 ) = 2 to 6 at incre- 
ments of 1. As a starting point, we computed 
a grid for an instantaneous burst with an IMF 
having a power-law slope of —2.35, M up = 100 
Mq, solar metallicity in stars and gas, and a sin- 
gle plane-parallel slab of gas with no dust; we will 
refer to this as model grid A. The solar abundance 
set was taken from Grevesse & Anders (1989) and 
Grevesse & Noels (1993). Other grids were com- 
puted as variations on this basic parameter set, 
with the following modifications made in different 
model runs: a constant star-formation rate; metal- 



licity 0.2, 0.4, or 2Zq in both stars and gas; and 
spherical geometry for the nebula. To assess the 
effects of the highest-mass stars, we also ran cus- 
tom model grids, via the STARBURST99 web 
site, with M up = 70 and 120 M Q . 

The depletion of heavy elements onto grains 
can result in marked changes to the emergent 
emission-line spectrum of an H II region, both 
by removing gas-phase coolants from the neb- 
ula and by grain absorption of ionizing photons, 
which will modify the effective shape of the ioniz- 
ing continuum. In metal-rich H II regions, these 
effects will tend to boost the strengths of the 
low-ionization emission lines relative to the dust- 
free case (Shields & Kennicutt 1995). To as- 
sess the effects of dust in transition nuclei, we 
calculated additional model grids which included 
dust grains with a Galactic ISM dust-to-gas ra- 
tio along with the corresponding gas-phase deple- 
tions. The dusty models were all calculated us- 
ing the solar abundance set for the undepletcd 
gas. Dust grains were assumed to have the optical 
properties of Galactic ISM grains, as described by 
Mathis, Rumpl, & Nordsieck (1977), Draine & Lee 
(1984), and Martin & Rouleau (1991). From the 
CLOUDY output files, we tabulated the strengths 
relative to FL3 of the major emission lines which 
are prominent in LINERs. 

The calculations were performed under the 
assumption that the H II region is ionization 
bounded. For this case, the outer extension of 
the cloud was set to be the radius at which T e 
falls to 4000 K, beyond which essentially no emis- 
sion is generated in the optical or UV lines. As 
a test, we ran a grid of models with the stopping 
temperature set to 1000 K, and we verified that 
the emission-line ratios were essentially identical 
to the default case of 4000 K stopping tempera- 
ture. We also verified that the important diag- 
nostic line ratios differed by < 0.1 dex between 
the spherical and plane-parallel cases when all 
other input parameters were unmodified, and all 
results discussed in this paper refer to the plane- 
parallel models. In the calculations, the longest 
timescales for atomic species to reach equilibrium 
were of order 10 3 years, much shorter than the 
evolution timescale of the stellar cluster, justi- 
fying the assumption that each time step of the 
cluster evolution could be used independently to 
calculate the nebular conditions. Table 1 gives a 



6 



summary of the model parameters, for the model 
grids which appear in the following discussion. 

4. Discussion 

4.1. Model Results 

The model results are displayed in Figures 1- 
9. To compare the model outputs with the ob- 
served properties of a variety of galaxy types, we 
have used the emission-line data compiled by Ho 
et al. (1997a). This catalog has the advantages of 
a homogeneous classification system, small mea- 
surement aperture (2" x 4"), and careful starlight 
subtraction to ensure accurate emission-line data. 
In order to reduce confusion and to keep the sam- 
ple of comparison objects to a reasonable number, 
we included only objects with unambiguous clas- 
sifications as H II, LINER, transition, or Seyfert. 
Objects with borderline or ambiguous classifica- 
tions, such as "LINER/ Seyfert," were excluded 
for clarity. The comparison sample was further 
reduced by excluding galaxies in which any of the 
emission lines Ha, H/3, [O III] A5007, [O I] A6300, 
[N II] A6583, or [S II] AA6716,6731 was undetected 
or was flagged as having a large uncertainty in flux 
("b" or "c" quality flags). The measured line ra- 
tios are corrected for both Galactic and internal 
reddening. 

Figure 1 plots the ratio [O III] A5007/H/3 
against [O I] A6300/Ha at a burst age of 4 Myr, 
for the solar-metallicity model grids A, B, C, and 
D. The model results are plotted for a density of 
njj = 10 3 cm~ 3 , as an approximate match to the 
density of n e = 600 cm -3 measured for NGC 4569 
(Ho et al. 1997a). The diagram shows that the 
instantaneous burst models (A and B) are a good 
match to the line ratios of the transition nuclei, for 
log U ~ —3.5. The dusty models from grid B fall 
more centrally within the region defined to con- 
tain transition objects, but the models without 
dust still closely match transition nuclei having 
lower [O I] A6300/Ha ratios. Figures 2 and 3 show 
the corresponding diagrams for [N II] A6583/Ha 
and for [S II] AA6716, 6731/Ha, respectively. In 
both cases we find that the single-burst models 
span the region occupied by transition nuclei in 
the diagnostic diagrams. 

In the constant star-formation rate models (C 
and D), the UV continuum remains softer than 
in the W-R-dominated phase of the instantaneous 



burst models, because of the ongoing formation of 
luminous O stars. As a result, the low-ionization 
emission lines are significantly weaker than in the 
instantaneous burst models at 4 Myr. These con- 
stant star-formation rate sequences are a reason- 
able match to the region of H II nuclei in the dia- 
gram, and for [N II] /Ha and [S II] /Ha the agree- 
ment with H II nuclei is improved at lower den- 
sities of n e = 10 2 cm -3 , a value more typical of 
H II nuclei. These constant star-formation rate 
models are probably appropriate for galaxies hav- 
ing spatially extended, ongoing star formation in 
their nuclei. 

We note that the model curves shown in the fig- 
ures should not be expected to follow the locus of 
H II nuclei in each plot, despite the fact that the 
models are generated with a starburst continuum. 
The range of line ratios observed in H II regions 
is primarily a sequence in metal abundance (Mc- 
Call, Rybski, & Shields 1985), while our models 
are shown as sequences in U for a given metallic- 
ity and density. Another point to note about the 
diagrams is that some of the transition galaxies 
fall outside the region nominally defined for tran- 
sition objects, particularly in the [O I] /Ha ratio 
(Figure 1). These galaxies were classified as tran- 
sition objects by Ho et al. (1997a) on the basis of 
meeting the majority of the classification criteria. 
Similarly, some overlap can be seen in the dia- 
grams between the regions occupied by LINERs 
and Seyfert nuclei; this again reflects the fact that 
galaxies span a continuous range in the values of 
these emission-line ratios. 

The variation of [O I] line strength as a func- 
tion of density is shown in Figure 4 for model grid 
A at t — 4 Myr. For the range of densities consid- 
ered (riH = 10 2 to 10 6 cm -3 ), the models closely 
overlap the transition region in the diagram at 
log U « —3.5. Introducing ISM depletion and 
dust grains to the nebula primarily increases the 
[O I] /Ha ratio at low density. As expected, the 
[O I] /Ha ratio increases with uh up to densities 
of 10 5 cm -3 , while at densities approaching the 
critical density of the A6300 transition (1.6 x 10 6 
cm~ 3 ) this ratio saturates and begins to turn over. 

Figure 5 shows the [O I] /Ha ratio as a function 
of burst age (at n e — 10 3 cm~ 3 ) for model grid A, 
for ages of 2 to 6 Myr. This diagram highlights 
the dramatic changes that W-R stars generate in 
the surrounding nebula. From 3 to 5 Myr after 
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the burst, when W-R stars dominate the UV con- 
tinuum, the harder ionizing continuum boosts the 
strength of [O I] by an order of magnitude and the 
model sequences appear adjacent to the transition 
region, with relatively little evolution in the line 
ratios during this period. As the burst ages be- 
yond 6 Myr, the [O I] /Ha ratio continues to fall, 
and the emission-line strengths drop rapidly as the 
ionizing continuum softens and its luminosity de- 
creases. At Z = 2Zq, the WR-dominated phase 
occurs slightly later, during time steps 4, 5, and 6 
Myr. The [N II] /Ha and [S II] /Ha ratios have a 
similar dependence on burst age, and are displayed 
in Figures 6 and 7, respectively; these results are 
quite similar to the calculations presented by Lei- 
therer ct al. (1992) to illustrate the effects of the 
W-R continuum on the [N II] /Ha ratio. While 
the [O I] /Ha ratios in the models at logC/ = —3.5 
are too low by ~ 0.1 — 0.2 dcx to fit within the 
nominal transition region, they still closely match 
those transition nuclei having relatively low values 
for [O I] /Ha, and the low-ionization line strengths 
can be further enhanced by the inclusion of dust 
and depletion (as in Figure 1). 

One puzzling aspect of Figure 5 is that for ages 
outside the range 3-5 Myr, the models predict 
[O I] /Ha ratios too low to match the majority 
of H II nuclei. Stasihska & Lcitherer (1996) and 
Martin (1997) discuss this same problem in the 
context of low-metallicity star burst galaxies. They 
propose that shocks generated by supernovae and 
stellar winds provide the additional [O I] emission, 
without making a significant contribution to the 
[O II] or [O III] line strengths. Shocks could play 
a similar role in transition galaxies, as in the model 
of Alonso-Herrero et al. (1999), although the lack 
of high-excitation UV line emission in transition 
nuclei is problematic for the shock hypothesis. A 
higher upper mass cutoff alleviates this problem 
to some extent, at least for very young bursts. In- 
creasing M up to 120 M Q boosts [O I] /Ha by - 0.2 
dex for ages of < 3 Myr. Due to the very short 
lifetimes of the highest-mass stars, however, the 
M up = 120, 100, and 70 M model grids result 
in identical emission-line spectra from ~ 4 Myr 
onward. 

Metal abundance is an additional parameter 
which must be considered. The models displayed 
up to this point were all calculated for a solar 
abundance set, while the nuclei of early- type spi- 



rals are likely to have enhanced heavy-element 
abundances. As discussed by Leitherer et al. 
(1999), the continuum shortwards of 228 A is 
strongest in the high-metallicity models, because 
the increased mass-loss rates lead preferentially to 
the formation of W-R stars at high metal abun- 
dance. Figure 8 plots the [O I] /Ha ratio for abun- 
dances of Z = 0.2, 0.4, 1, and 2 Z Q , at a den- 
sity of jih = 10 3 cm~ 3 and t = 4 Myr. From 
this diagram, it is clear that solar or higher abun- 
dances are necessary to match the [O I] strengths 
of the transition nuclei; at lower abundances the 
line ratios are a better match to those of the high- 
excitation (low-metallicity) H II nuclei. Figure 9 
displays the density dependence of the [O I] /Ha 
ratio for the Z = 2Zq model grid; by comparison 
with the solar-metallicity model grid in Figure 4, 
the higher abundances result in a lower-excitation 
spectrum with enhanced [O I] emission, due to the 
harder extreme-UV continuum. 

It would be advantageous to compare the model 
results with a wider variety of emission lines. Un- 
fortunately, measurements of other optical emis- 
sion lines are scarce for transition nuclei. The Ho 
ct al. survey did not include the [O II] A3727 line, 
and there is no other homogeneous catalog of [O II] 
measurements for transition galaxies. To be con- 
sistent with a LINER or transition-type classifica- 
tion, a model calculation must result in the flux 
ratio [O II] A3727 / [O III] A5007 > 1. In fact, all 
of the models with log U < — 3 do satisfy this crite- 
rion. Thus, any of our models which is consistent 
with the Ho et al. LINER or transition classifi- 
cation criteria is also consistent with the original 
Heckman (1980) criterion for the [O II] /[O III] ra- 
tio in LINERs. The relative strengths of UV lines 
such as C II] A2326, C III] A1909, and C IV A1549 
can provide further diagnostics, but none of these 
lines is detected in NGC 4569 (Maoz et al. 1998). 
The only other transition nucleus having HST UV 
spectra available is NGC 5055, and its spectrum 
appears to be devoid of UV emission lines as well 
(Maoz et al. 1998). 

We ran one additional model grid to test 
whether different model atmospheres for O stars 
would lead to different results. The STAR- 
BURST99 continua were calculated using stel- 
lar atmosphere models compiled by Lejeune et 
al. (1997), which are based on the Kurucz (1992) 
model set for the massive stellar component. The 
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recent CoStar model grid of Schaerer & de Koter 

(1997) , which includes non-LTE effects, stellar 
winds, and line blanketing for O stars, makes 
dramatically different predictions for the ioniz- 
ing spectra. As shown by Schaerer & Vacca 

(1998) , the CoStar models yield a luminosity in 
the He ++ continuum which is four orders of mag- 
nitude greater than that predicted by the Kurucz 
models, for the most massive O stars which domi- 
nate the UV luminosity at burst ages of < 3 Myr. 
In the CoStar-based models the photon output 
of the cluster below 228 A is essentially constant 
from to 5 Myr. To investigate the effects of 
this harder O-star continuum on the emission-line 
spectra, we ran a grid of models using the evolving 
starburst continua computed by Schaerer & Vacca 
(1998) with the CoStar atmospheres. Model pa- 
rameters were the same as for model grid A except 
that an upper mass limit of 120 Mq was used. We 
find that using the CoStar atmospheres has a rel- 
atively minor effect on our results. In comparison 
with the STARBURST99-based models having 
M up = 120M Q , the CoStar model grid yields an 
increase in the [O I] /Ha ratio of <~ 0.1 — 0.15 dex 
for t < 6 Myr, while [N II] /Ha and [S II] /Ha are 
essentially unaffected. During the period t < 3 
Myr, the CoStar-based models result in an H II 
region spectrum, demonstrating that W-R stars 
are still required in order to generate LINER or 
transition- type line ratios. 

The strength of the [Ca II] emission lines at 
7291 and 7324 A is often used as a diagnostic 
of dust and depletion, because in the absence of 
depletion these lines are predicted to be strong 
in photoionized gas (e.g., Kingdon et al. 1995; 
Villar-Martm & Binette 1996). (The A7291 line 
is a cleaner diagnostic since A7324 is blended with 
[O II] A7325.) However, for a 4 Myr-old burst 
with nebular conditions of n e = 10 3 cm~ 3 , log U = 
—3.5, and an undcplcted solar abundance set, our 
calculations yield a maximum prediction of only 
0.2 for the ratio of [Ca II] AA7291, 7324 to H/3. 
Only at very low ionization parameters (< 10 -4 5 ) 
does the [Ca II] emission become stronger than 
H/3. Since H/3 is only barely visible in the spec- 
tra of many transition objects (prior to careful 
starlight subtraction, at least), typical observa- 
tions may not have sufficient sensitivity to detect 
faint [Ca II] lines in these objects. High-quality 
spectra of LINERs do not show [Ca II] emission 



(Ho et al. 1993), indicating that Ca is likely to be 
depleted onto grains in these objects, but similar 
data are not generally available for transition nu- 
clei. If the [Ca II] lines are found to indicate a 
high level of depletion onto dust grains in transi- 
tion nuclei, this would also provide a further argu- 
ment against shock-heating models, as shocks will 
tend to destroy grains (e.g., Morse et al. 1996). 

4.2. The Nature of Transition Nuclei 

The results shown in the preceding figures 
demonstrate that the starburst models are in fact 
able to reproduce the major diagnostic emission- 
line ratios of transition nuclei with reasonable ac- 
curacy, during the period t = 3-5 Myr when W-R 
stars are present. For a density of 10 3 cm~ 3 and 
an age of 4 Myr, the solar-metallicity models with 
and without depletion bracket the range of values 
observed in real transition nuclei for the line ra- 
tios [O I]/Ha, [N II]/Ha, and [S II]/Ha. We do 
not attempt to fine-tune a model to produce an 
exact match with the spectrum of NGC 4569, but 
the basic solar-metallicity dust-free model at t = 4 
Myr with nn = 10 3 cm~ 3 and log U = —3.5 closely 
fits the observed [O I] /Ha ratio, while overpredict- 
ing [S II] /Ha by ~ 0.2 dex and underpredicting 
[N II] /Ha by - 0.1 dex. 

We emphasize that the starburst models are 
only able to produce transition-type spectra for 
the case of an instantaneous burst; that is, when 
the burst duration is shorter than the timescale 
for evolution of the most massive stars. Multiple- 
burst populations can only yield a transition spec- 
trum if the dominant population is ~ 3 — 5 Myr 
old and the older or younger bursts do not con- 
tribute significantly to the ionizing photon bud- 
get. Models with a constant star- formation rate 
produce H II region spectra at all ages, as the 
softer ionizing continua do not produce sufficient 
[O I] A6300 emission in the surrounding H II region 
to match transition-type spectra. The parameter 
which is most important for determining the hard- 
ness of the ionizing continuum is the number ratio 
of W-R stars to O stars, which exceeds ~ 0.15 
during the W-R-dominated phase in the STAR- 
BURST99 models at solar metallicity, and ap- 
proaches or exceeds unity at Z = 2Zq. In the con- 
stant star-formation rate models at solar metal- 
licity, the W-R/O ratio levels off at <~ 0.06 after 
about 4 Myr. The compact size of the NGC 4569 
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nucleus is consistent with the requirement that the 
burst duration must be brief (< 1 Myr) in order to 
generate a transition-type spectrum. The FWHM 
size of the starburst core in NGC 4569 is only ~ 10 
pc. For such a burst to occur in < 1 Myr would 
require a propagation speed for star formation of 
only ~ 10 km s _1 . In fact, the typical velocities in 
the NGC 4569 nucleus are much greater than 10 
km s _1 : the [N II] A6583 line has a velocity width 
of 340 km s- 1 (Ho ct al. 1997a). Thus, the NGC 
4569 nucleus could represent the result of a single, 
rapid burst of star formation. 

Although our results suggest that transition 
galaxy spectra may be attributed to a starburst 
with a high W-R/O-star ratio, the demographics 
of transition nuclei and H II nuclei indicate that 
many transition galaxies are probably not formed 
by this mechanism. In the STARBURST99 mod- 
els, the W-R-dominated phase in an instantaneous 
burst lasts for ~ 3 Myr (i.e., 3 time steps in the 
calculations). An H II region surrounding an in- 
stantaneous burst will be visible for ~ 6 Myr, af- 
ter which the emission lines will fade rapidly (e.g., 
Garcia- Vargas & Diaz 1994). Thus, for an in- 
stantaneous burst population, the transition phase 
and the H II nucleus phase will have approxi- 
mately equal lifetimes. If all H II nuclei con- 
sisted of instantaneous burst stellar populations 
with nebular conditions conducive to the forma- 
tion of transition-type spectra, then H II nuclei 
and transition nuclei should be roughly equal in 
number. In reality, it is likely that a large fraction 
of star-forming nuclei contain multiple bursts of 
star formation and/or conditions of low density or 
low metallicity, so all star-forming nuclei should 
not be expected to evolve through a transition- 
type phase. Although it is difficult to make spe- 
cific predictions, it is probably safe to conclude 
that for a given Hubble type, transition nuclei gen- 
erated solely by starbursts should be considerably 
less numerous than ordinary H II nuclei. 

The statistics compiled by Ho et al. (1997b) 
provide a basis for comparison. In early-type 
galaxies (E and SO), transition nuclei outnumber 
H II nuclei by a 3-to-l margin. Only for Hubble 
types Sb and later do H II nuclei begin to outnum- 
ber transition nuclei by a factor of 2 or more. The 
most straightforward interpretation of this trend 
is that in early-type host galaxies, the majority 
of transition nuclei are actually AGN/H II region 



composites, as proposed by Ho et al. (1993) and 
others. At intermediate and late Hubble types, the 
population of transition nuclei may consist of both 
composite objects and "pure" starbursts evolving 
through the W-R-dominatcd phase. 

The presence of transition nuclei in a small frac- 
tion (~ 10%) of elliptical galaxies (Ho et al. 1997b) 
presents a particularly intriguing problem. The 
Ho ct al. survey detected five transition nuclei 
in ellipticals but not a single case of an ellipti- 
cal galaxy hosting an H II nucleus. Given that 
the models which have been considered for transi- 
tion nuclei involve star formation, cither alone or 
in combination with an AGN, this observation is 
rather puzzling. Perhaps faint AGNs in elliptical 
nuclei can produce transition-type spectra with- 
out substantial star formation activity. Four of 
the five transition nuclei found in ellipticals by Ho 
et al. (1997a) have borderline or ambiguous spec- 
troscopic classifications, however, so "pure" tran- 
sition nuclei in ellipticals are evidently quite rare. 

Given these results, one might expect to see 
transition-type emission spectra in some fraction 
of disk H II regions in spiral galaxies, but in fact 
such spectra are never found. Single-burst models 
for disk H II region spectra are only compatible 
with observed line ratios for model ages of t < 3 
Myr (Brcsolin et al. 1999), as the harder ioniz- 
ing spectrum after 3 Myr makes the models over- 
predict the strengths of the low-ionization lines. 
Bresolin et al. suggest that either current stellar 
evolution models are at fault, or that disk H II 
regions are disrupted before reaching an age of 3 
Myr, in which case the W-R phase would not be 
observed in the nebular gas. An alternate (and 
perhaps more attractive) possibility is that the 
majority of H II nuclei, as well as disk H II regions, 
are better described by the models with constant 
star-formation rate, or contain multiple bursts of 
star formation with an age spread of a few Myr, 
which would result in a spectrum similar to the 
constant star-formation rate models. 

For understanding the physical nature of tran- 
sition objects, the observational challenge is to 
search for any unambiguous signs of nonstellar ac- 
tivity. Detection of broad Ha emission, or a com- 
pact source of hard X-ray emission with a power- 
law spectrum, would provide evidence for an AGN 
component. High- resolution optical spectra (from 
HST) could provide a means to spatially resolve 
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a central AGN-dominated narrow-line region from 
the surrounding starburst-dominated component. 
Since direct evidence for accretion-powered nu- 
clear activity in transition nuclei is generally lack- 
ing, it should not be assumed that any given tran- 
sition object actually contains an AGN unless ob- 
servations specifically support that interpretation. 

One further effect that should be considered 
in starburst models in the future is photoioniza- 
tion by the X-rays generated by the starburst. X- 
ray binaries and supernova remnants will provide 
high-energy ionizing photons, resulting in a spa- 
tially extended source of soft X-ray emission as 
observed in the nucleus of NGC 4569 (Terashima 
et al. 1999), for example. (The massive main 
sequence stars will contribute only a negligible 
amount to the total X-ray luminosity of a star- 
burst; see Helfand & Moran 1999.) Photoioniza- 
tion by X-rays will naturally lead to an enhance- 
ment of the low-ionization forbidden lines, and this 
could contribute to the excitation of some transi- 
tion galaxies. 

4.3. LINERs 

The strength of [O I] A6300 is the key distin- 
guishing factor between LINERs and transition 
nuclei, and matching the observed strength of this 
line is the major challenge of starburst models 
for LINERs. Our calculations show that LINER 
spectra can only be generated by the STAR- 
BURST99 clusters under a very specific and lim- 
ited range of circumstances. Model grids A and 
B, while matching the [O I] / Ha ratio of transi- 
tion nuclei quite well, do not overlap at all with 
the main cluster of LINERs in Figure 1, even at 
high densities and even when depletion and dust 
grains arc included. Only grid G with Z = 2Z & 
is able to replicate the high [O I] / Ha ratios of 
most LINERs, and only during t w 4-6 Myr and 
at densities of «h ^ 10 5 cm~ 3 . In agreement with 
previous models, we find that values of log U « 
—3.5 to —3.8 reproduce the observed [O I] / Ha 
ratios of LINERs. However, at such high densities 
the models underprcdict the strengths of [S II] and 
[N II] relative to Ha. Single-zone models require 
n e < 10 5 cm~ 3 to match the [N II] /Ha ratios of 
LINERs and n e < 10 4 cur 3 for [S II] /Ha. 

Agreement with LINER spectra can be achieved 
with a simple two-zone model, in which high- 
density and low-density components are present, 



similar to the scenario proposed by Shields (1992). 
As an example, a two-component model con- 
structed from grid A containing gas at (n e = 10 3 
cm~ 3 , U = 10~ 35 ) and at (n e = 10 5 cm~ 3 , 
U = 10~ 4 ) produces emission-line ratios which 
are consistent with all the LINER classification 
criteria of both the Heckman (1980) and Ho et al. 
(1997a) systems, if the two density componenets 
are scaled so as to contribute equally to the total 
H/3 luminosity. As a local comparison, obser- 
vations of near-infrared Fe II emission indicate 
the presence of clouds having n e > 10 5 cm~ 3 in 
the Galactic center region (DePoy 1992), so it is 
plausible that other galactic nuclei may contain 
ionized gas at similarly high densities even in the 
absence of an observable AGN. A starburst origin 
for some LINER 2 nuclei would provide a natural 
explanation for the lower values of the X-ray/Ha 
flux ratio seen in these objects, in comparison with 
AGN-like LINER 1 nuclei (Terashima et al. 1999). 

It seems unlikely, however, that many LINERs 
arc generated by this starburst mechanism. About 
15% of LINERs are known to have a broad com- 
ponent of the Ha emission line, indicating a prob- 
able AGN (Ho et al. 1997b). By analogy with 
the Seyfert population, a much larger fraction of 
LINERs is likely to have broad-line regions which 
are either obscured along our line of sight, or are 
simply too faint to be detected in ground-based 
spectra against a bright background of starlight. 
Many LINERs show signs of nuclear activity that 
cannot be explained by stellar processes: compact 
flat-spectrum radio sources or jets, compact X-ray 
sources with hard power-law spectra, or double- 
peaked broad Balmer-line emission, for example. 
As an increasing body of observational work sup- 
ports the idea that many LINERs are in fact 
AGNs, there is less incentive to consider purely 
stellar models for their excitation. 

Demographic arguments, similar to those given 
above for transition nuclei, can be applied for the 
LINER population. Since the LINER phase only 
occurs for instantaneous bursts at high density 
and high metallicity, the starburst scenario im- 
plies that H II nuclei should be considerably more 
numerous than starburst-generated LINERs for a 
given Hubble type. While LINERs are common 
in early- type hosts, H II nuclei are not found in 
elliptical hosts and are seen in fewer than 10% 
of SO galaxies (Ho et al. 1997b). This dispar- 
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ity is a strong argument against a starburst ori- 
gin for those LINERs in early-type galaxies. In 
later Hubble types the situation is less clear, how- 
ever. H II nuclei occur in ~ 80% of spirals of type 
Sc and later, while LINERs occur in just 5% of 
these galaxy types (Ho et al. 1997b). It is con- 
ceivable that some of the LINERs in intermediate 
to late-type hosts could have a starburst origin, 
and this issue could be resolved by further UV 
and X-ray observations in the future. Interest- 
ingly, the Ho et al. survey did not find any exam- 
ples of broad Ha emission in LINERs or transi- 
tion objects with hosts of type Sc or later; perhaps 
star formation plays a more prominent role than 
accretion-powered activity in these objects. 

While a few LINERs show spectral features of 
young stars in the UV (Maoz et al. 1998), the qual- 
ity of the observational data is poor in comparison 
with the NGC 4569 UV spectrum, and it is diffi- 
cult to set meaningful constraints on the age of 
the young stellar population. NGC 404 is a pos- 
sible candidate for a starburst-generated LINER, 
but in its UV spectrum the P Cygni features are 
weak in comparison with NGC 4569, indicating ei- 
ther an older burst population or dilution by a fea- 
tureless AGN continuum (Maoz et al. 1998). The 
LINERs having UV spectral features from massive 
stars may also host obscured AGNs which can be 
detected in other wavebands. For example, the 
UV continuum of the LINER NGC 6500 appears 
to have its origin in hot stars (Barth et al. 1997; 
Maoz et al. 1998), but observations of a parsec- 
scale radio jet unambiguously demonstrate that 
nonstellar activity is occurring as well (Falcke et 
al. 1998). 

4.4. W-R Galaxies with LINER or Transition- 
Type Spectra 

The starburst models presented here run into 
two obvious problems. First, W-R galaxies are 
almost never known to have LINER or transition- 
type spectra. Second, LINERs and transition nu- 
clei almost never show W-R features in their spec- 
tra. Is there any way to reconcile the starburst 
models with these facts? 

W-R galaxies are identified by the appearance 
of the 4650 A blend in their spectra (e.g., Kunth & 
Sargent 1981). Since the formation of W-R stars 
is enhanced at high Z, the strength of this feature 
relative to H/3 increases dramatically with metal- 



licity, from < 0.1 at Z < 0.4Z© to - 0.5 - 4 at 
Z > Z Q (Schaerer & Vacca 1998). However, in 
the nuclei of early-type spirals where high metal- 
licities are expected, a nuclear starburst will be 
surrounded by the old stellar population of the 
galactic bulge, making the detection of the W-R 
bump extremely difficult (Mas-Hesse et al. 1999). 
Most of the currently known W-R galaxies are 
late-type spirals or irregular galaxies (Schaerer et 
al. 1999) in which the W-R bump is visible against 
the nearly featureless starburst continuum. When 
the W-R bump is detected in H II galaxies, its am- 
plitude above the continuum level is generally far 
smaller than that of H/3 or even H7 (e.g., Kunth 
& Sargent 1981). In most of the LINER and tran- 
sition galaxy spectra in the catalog of Ho et al. 
(1995), however, H/3 barely appears and H7 is 
too weak to be visible at all prior to continuum 
subtraction. Even in a high-mctallicity environ- 
ment, where the total intensity of the W-R bump 
can be comparable to that of H/3, the amplitude 
of the W-R bump above the continuum will be 
much lower than that of H/3 because the flux in the 
W-R feature is spread over ~ 70 - 100 A. Thus, 
the detection of W-R emission in galactic nuclei is 
strongly biased toward late-type, bulgcless galax- 
ies. In late- type or dwarf irregular galaxies where 
the W-R bump is visible, the W-R/O-star ratio is 
expected to be much smaller owing to the lower 
metallicity, and the resulting softer ionizing spec- 
trum will tend to produce an H II region spectrum 
rather than a transition object. The gas density 
as a function of Hubble type may play a role as 
well; in a study of H II nuclei, Ho et al. (1997c) 
find a weak trend toward lower nebular densities 
in later- type host galaxies. 

Observational detection of W-R features in 
transition nuclei is perhaps the clearest test of 
the starburst models, if sufficiently sensitive ob- 
servations can be obtained. The UV spectrum of 
NGC 4569 is consistent with an age of ~ 3 — 6 
Myr, an age at which W-R stars are expected 
to be present. Previous optical spectra have not 
revealed the 4650 A W-R bump in NGC 4569, 
but further observations with high S/N and small 
apertures would be worthwhile. The lack of He II 
A1640 emission in the UV spectrum of NGC 4569 
is potentially a more serious problem since the 
burst population should dominate at short wave- 
lengths. The models of Schaerer & Vacca (1998) 
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predict an equivalent width of at least 2 A in 
the W-R-generated A1640 line during the period 
3-6 Myr for an instantaneous burst of solar or 
higher mctallicity, while the observed upper limit 
of /(1640) < 2.0 x lCT 15 erg s" 1 craT 2 (Maoz et 
al. 1998) corresponds to an equivalent width limit 
of < 0.3 A. It should be noted, however, that the 
A 1640 line lies at the extremely noisy blue end 
of the FOS G190H grating setting, in a region 
where detection of emission or absorption features 
is difficult. 

Two W-R galaxies may provide useful points of 
reference for the starburst models. NGC 3367 is 
classified by Ho et al. (1997a) as an H II nucleus 
on the basis of its [O I] /Ha and [S II] /Ha ratios, 
although its [N II] /Ha ratio of 0.83 is more con- 
sistent with a LINER or transition-type classifica- 
tion and its emission lines are markedly broader 
than those of typical H II nuclei. Alonso-Herrero 
et al. (1999) describe NGC 3367 as a starburst- 
dominated transition object (see also Dekker et al. 
1988). The 4650 A W-R bump was noted by Ho 
et al. (1995), who also suggested a LINER/H II 
classification and a composite source of ioniza- 
tion. As a borderline H II nucleus/transition ob- 
ject with clear evidence for W-R stars, this object 
deserves further study, to determine whether there 
is indeed an AGN or whether the enhanced low- 
ionization emission may be the result of ionization 
by the W-R population. The electron density of 
835 cm -3 measured from the [S II] doublet (Ho et 
al. 1997a) is also noteworthy, as this is among the 
highest densities found for an H II nucleus in the 
Ho ct al. survey. 

Another intriguing object is the nucleus of NGC 
6764, which has been classified variously as a 
Seyfert, a LINER, and a starburst by different 
authors (see Goncalves et al. 1999). This galaxy 
exhibits prominent emission in the 4650 A W-R 
blend (Osterbrock & Cohen 1982). A recent study 
by Eckart et al. (1996) demonstrates that the nar- 
row emission lines are consistent with a LINER 
classification, but there are no unambiguous signs 
of nonstellar activity in the nucleus. Eckart et al. 
(1996) find that the nucleus contains ~ 3600 W-R 
stars, and that the overall properties of the ob- 
ject are consistent with ionization by the starburst 
alone, rather than by a starburst/ AGN composite. 
If this conclusion is confirmed by further obser- 
vations, NGC 6764 could be considered the best 



candidate for a LINER photoionized by a star- 
burst during its W-R-dominated phase. Alonso- 
Herrero et al. (1999) derived an age of 9-10 Myr 
for the starburst in NGC 6764 based on near- 
infrared emission-line diagnostics, but this age is 
inconsistent with the presence of W-R stars, at 
least for the case of an instantaneous burst. 

Compared with typical LINERs and transition 
nuclei, conditions for detection of W-R spectral 
features in these two objects are perhaps more fa- 
vorable. Both host galaxies are of late Hubble 
types (SBc for NGC 3367 and SBbc for NGC 6764) 
in comparison with the majority of LINERs and 
transition nuclei, so the level of contamination by 
the surrounding old stellar population is relatively 
low. Furthermore (and partly as a result of this), 
the emission-line equivalent widths in these two 
nuclei are relatively high for LINERs or transition 
nuclei. If the emission-line spectra of NGC 3367 
or NGC 6764 had been superposed on a luminous 
early-type spiral bulge, the W-R emission might 
never have been noticed. 

The detection of W-R emission in a LINER 
does not automatically imply a purely starburst 
origin for the emission lines, of course, since star- 
bursts and AGNs are often known to coexist. Mrk 
477 is a well-known example of a Seyfert 2 galaxy 
having a large population of W-R stars in its nu- 
cleus (Hcckman et al. 1997). In the context of the 
starburst model, however, the crucial test is to 
search for additional examples of LINERs or tran- 
sition nuclei which exhibit a high ratio of W-R to 
O stars but no signs of accretion-powered activity. 

4.5. Caveats and Limitations 

The most important limitation of these calcu- 
lations comes from the accuracy of the input con- 
tinua. The conclusion that the starburst models 
are able to reproduce transition or LINER spec- 
tra under some circumstances depends crucially on 
the presence of W-R stars to provide a hard and 
luminous ionizing continuum. Unfortunately, the 
continuum shape and luminosity of W-R stars in 
the extreme-UV band are quite uncertain, partic- 
ularly in the Hc ++ continuum where stellar winds 
have a dramatic effect. Several up-to-date reviews 
of the numerous difficulties involved in modeling 
W-R spectra can be found in the volume edited 
by van der Hucht et al. (1999). The STAR- 
BURST99 model grid uses the W-R atmospheres 
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of Schmutz ct al. (1992), which are calculated for 
a pure helium composition, but more recent at- 
mosphere models are beginning to include the ef- 
fects of line-blanketing, as well as clumping and 
departures from spherical symmetry. As discussed 
by Leitherer et al. (1999), the W-R/O-star ra- 
tio is also extremely model dependent, and may 
be revised in future generations of models. This 
would have a direct effect on the strength of the 
cxtrcmc-UV continuum and consequences for the 
nebular emission lines. Furthermore, the STAR- 
BURST99 models neglect binary evolution, al- 
though this is more likely to affect the W-R/O 
ratio at low metallicity. Since the results of the 
transition-object models are highly dependent on 
the most uncertain portion of the W-R spectrum, 
new photoionization calculations should be com- 
puted to assess the impact of different W-R evo- 
lution and atmosphere models in the future. 

The shape of the IMF in starburst regions is 
a subject of some debate, and the possible vari- 
ation of the IMF with metallicity is of particular 
importance for galactic nuclei, which are likely to 
have Z > Z Q . Kahn (1974) and Shields & Tins- 
ley (1976) suggested that M up should be lower in 
regions of higher metal abundance, but this issue 
has not been settled definitively. Star-count obser- 
vations demonstrate that the IMF slope and M up 
do not appear to vary with metallicity (Massey 
et al. 1995), at least for Z < Zq. While nebular 
diagnostics in H II galaxies are generally consis- 
tent with a Salpeter IMF with M up w 100 M Q 
at subsolar metallicity (e.g., Stasihska & Leitherer 
1996), at high metallicity the observational sit- 
uation is somewhat ambiguous. Bresolin ct al. 
(1999) find that the mean stellar temperature in 
H II regions decreases significantly with increas- 
ing Z, and that the He I A5876 / H/3 ratios of 
H II regions at Z « 2Z Q are more consistent 
with M up = 30 M Q than with M up = 100 M . 
Such a low value for M up would pose serious dif- 
ficulties for any starburst models of LINERs and 
transition nuclei, as the massive progenitors of W- 
R stars would not be present. Counterbalancing 
this trend, the strong tidal forces, turbulence, and 
magnetic field strengths in galactic nuclei may act 
to raise the Jeans mass and favor the formation of 
more massive stars (Morris 1993). In the Galac- 
tic center, there are stars with initial masses of 
~ 100 M Q (Krabbe et al. 1995), and one Galactic 



center object (the Pistol star) may have Mi„j t i a i 
as high as 200-250 M Q (Figer et al. 1998). Thus, 
the proposed trend toward lower values of M up at 
high metallicity in disk H II regions may not ap- 
ply to galactic nuclei. Detailed comparison of the 
UV spectra of galaxies such as NGC 4569 with 
starburst population synthesis models can provide 
useful constraints on the population of high-mass 
stars in nuclear starbursts. 

Despite these uncertainties, these photoioniza- 
tion models have a major advantage compared 
with previous generations of W-R or O-star mod- 
els for LINERs and transition nuclei, in that the 
STARBURST99 models with standard parame- 
ters are constructed to represent the actual stellar 
populations in starbursts, to the best of current 
knowledge. Previous O-star models (Filippenko 
& Terlevich 1992; Shields 1992) required the pres- 
ence of hypothetical, unusually hot stars in order 
to explain LINER or transition spectra, and they 
did not address the evolution of the young stellar 
population at all. The starburst models presented 
here provide a more plausible mechanism to gener- 
ate a transition-type spectrum, even if this model 
may apply only to a relatively small fraction of the 
population of transition galaxies. 

5. Conclusions 

Our primary conclusion is that for standard 
starburst parameters and for nebular conditions 
which may be typical of galactic nuclei, the star- 
burst models are able to reproduce the important 
diagnostic emission-line ratios for LINER/H II 
transition galaxies, otherwise known as weak-[0 I] 
LINERs. The key ingredient needed to generate a 
transition-type spectrum is a UV continuum domi- 
nated by W-R stars, a condition which occurs dur- 
ing t = 3-5 Myr after an instantaneous burst. A 
transition-type emission spectrum may thus be a 
phase in the evolution of some nuclear H II re- 
gions in which the ionizing continuum is gener- 
ated by a single-burst stellar population. The 
models are also able to produce an [O I] / Hen 
ratio high enough to match LINER spectra, but 
only for conditions of above-solar metallicity com- 
bined with the presence of high-density (> 10 5 
cm~ 3 ) clouds. A sensitive search for W-R spec- 
tral features in transition nuclei would provide a 
test of this starburst scenario. This model may 
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apply only to a small fraction of LINERs and tran- 
sition nuclei; many LINERs and some transition 
objects show clear signs of nonstellar activity, and 
the starburst models may not apply at all to ob- 
jects in early- type host galaxies. Further multi- 
wavelength observations of transition nuclei will 
be of great utility for determining what fraction 
of them contain genuine active nuclei, and what 
fraction appear to be purely the result of stellar 
phenomena. 
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Table 1 
Model Parameters 



Model Grid SF Law 


z/z & 


Dusty? 
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I 
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C 
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I 
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N 
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I 


2.0 


N 



Note. — All models listed above are calculated for M up 
= 100 Mq, IMF power-law slope of —2.35, and for plane- 
parallel nebular geometry. The star-formation law in Col- 
umn 2 denotes I = instantaneous burst, and C = constant 
star-formation rate. The metallicity Z refers to the stars 
and to the undepleted abundance of the nebular gas. 
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Figure Captions 

Fig. 1.- 

Line-ratio diagram of [0 III] A5007 / H/3 against 
[O I] A6300 / Ha, for model grids A (solid line), 
B (long-dashed line), C (short-dashed line), and 
D (dot-dashed line), at a burst age of 4 Myr and 
n e = 10 3 cm~ 3 . The input continuum has solar 
metallicity, IMF power-law slope —2.35, and M up 
= 100 Mq. The following description applies to 
this and all subsequent plots: the small squares 
along each model line correspond to the model grid 
points at log U = —4, —3.5, —3, —2.5, and —2, with 
U increasing upward along the line. The points 
plotted represent galaxies from the Ho et al. cata- 
log, as follows: Small circles: H II nuclei. Squares: 
LINER/H II transition objects. NGC 4569 is rep- 
resented by an open square. Triangles: "pure" 
LINERs. Crosses: Seyfert nuclei. The dotted line 
encloses the region defined for transition nuclei ac- 
cording to the criteria of Ho et al. (1997a). 

Fig. 2. — Same as Figure 1, but for the line ra- 
tio [N II] A6583 / Ha. The dotted line encloses 
the region occupied by both LINERs and transi- 
tion objects, according to the criteria of Ho et al. 
(1997a). 

Fig. 3. — Same as Figure 1, but for the line ratio 
[S II] AA6716, 6731 / Ha. The dotted line encloses 
the region occupied by both LINERs and transi- 
tion objects, according to the criteria of Ho et al. 
(1997a). 

Fig. 4.— Density effects on the [O I] A6300 line 
strength, for model grid A. Plot symbols for in- 
dividual data points are the same as for Figure 1. 
The model sequences are shown for age 4 Myr and 
nebular density of log (nn/cm -3 ) = 2, 3, 4, 5, and 
6. 

Fig. 5.— Burst age effects on the [O I] A6300 
line strength, for model grid A. Plot symbols are 
the same as for Figure 1. The model sequences 
are shown for nebular density n e = 10 3 cm -3 and 
burst ages of 2, 3, 4, 5, and 6 Myr. 

Fig. 6.— Burst age effects on the [N II] A6583 
line strength, for model grid A. Plot symbols are 
the same as for Figure 2. The model sequences 
are shown for nebular density n e = 10 3 cm~ 3 and 



burst ages of 2, 3, 4, 5, and 6 Myr. 

Fig. 7.— Burst age effects on the [S II] AA6716, 
6731 line strength, for model grid A. Plot symbols 
are the same as for Figure 3. The model are shown 
for nebular density n e — 10 3 cm' 3 and burst ages 
of 2, 3, 4, 5, and 6 Myr. 

Fig. 8.— Metallicity effects on the [O I] A6300 / 
Ha flux ratio. Plot symbols are the same as for 
Figure 4. Model grids are plotted for the cases 
of Z/Zq = 0.2 (dot-dashed line; from grid E), 0.4 
(long-dashed line; from grid F), 1 (solid line; from 
grid A), and 2 (short-dashed line; from grid G). 
Each sequence is shown for «h = 10 3 cm~ 3 and 
t = 4 Myr. 

Fig. 9, The [O I] A6300 / Ha flux ratio for the 
case of Z = 2Z Q (model grid G), at t = 4 Myr. 
Plot symbols are the same as for Figure 4. The 
model sequences are plotted for densities of log 
(nj/cm -3 ) = 3, 4, 5, and 6. 
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